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We investigate theoretically the possibility of a wetting transition induced by geometric roughness 
of a solid substrate for the case where the flat substrate does not show a wetting layer. Our 
approach makes use of a novel closed-form expression which relates the interaction between two 
sinusoidally modulated interfaces to the interaction between two flat interfaces. Within the harmonic 
approximation, we find that roughness-induced wetting is indeed possible if the substrate roughness, 
quantified by the substrate surface area, exceeds a certain threshold. In addition, the molecular 
interactions between the substrate and the wetting substance have to satisfy several conditions. 
These results are expressed in terms of a lower bound on the wetting potential for a flat substrate 
in order for roughness-induced wetting to occur. This lower bound has the following properties: A 
minimum is present at zero or very small separation between the two interfaces, as characteristic 
for the non-wetting situation in the flat case. Most importantly, the wetting potential needs to have 
a pronounced maximum at a separation comparable to the amplitude of the substrate roughness. 
These findings are in agreement with the experimental observation of roughness-induced surface 
premelting at a glass-ice interface as well as the calculation of the dispersion interaction for the 
corresponding glass-water-ice system. 
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I. INTRODUCTION 

The phenomenon of wetting has been the subject of intense attention and a fairly good understanding of the basic 
concepts and mechanisms has emerged In the simplest case, a solid inert surface is put in contact with an 

under-saturated vapor of a second substance. Typically, a molecularly thick liquid-like film will form on the substrate 
surface due to favorable molecular interactions. The liquid film is in equilibrium with its under-saturated vapor, thus 
giving rise to a second interface (the emerging liquid- vapor interface), referred to hereafter as the liquid interface. 
Depending on the detailed molecular interactions between all three phases and the resulting interfacial energies, the 
liquid film can either grow to macroscopic thickness or remain finite as coexistence between the liquid and its vapor 
is approached. The first situation corresponds to complete wetting with a diverging film thickness, the second case is 
called incomplete or partial wetting. 

Two other, closely related situations are possible: (i) The liquid layer (e.g., water) on the inert substrate can be 
in equilibrium with its solid phase (ice) at temperatures below the melting point. In this case the vapor is replaced 
by a solid and the appearance of a thin liquid layer between the substrate and the solid phase indicates interfacial 
premelting. Note that the third phase (the ice) is entirely different from the solid substrate, (ii) The substrate itself 
can be a solid in equilibrium with its vapor phase. Here the formation of a thin liquid layer as three-phase coexistence 
is approached corresponds to surface premelting ^ . The phenomenological description of these scenarios does not 
differ from the wetting situation, and one finds the analogous phenomena corresponding to partial and complete 
wetting. 

In early theoretical studies, the solid substrate was assumed to be flat and homogeneous. However, in most 
experimental and technological situations the substrate is both rough and inhomogeneous. For complete wetting 
upon approaching coexistence, where the liquid forms a thin and continuous film, the influence of substrate roughness 
and chemical disorder has recently been investigated theoretically and experimentally [p~0| — |l3|| in great detail. 

It was found that heterogeneity and roughness of the solid substrate in conjunction with long-range van der Waals 
interactions cause equilibrium undulations of the liquid film surface. Surface tension, on the other hand, acts as a 
damping mechanism which reduces the amplitude of undulations for thicker films. The theoretical results (^^] were 
verified recently in small-angle x-ray scattering |ll]-[l3| . 

Yet another realization of the wetting phenomenon is obtained if a non-volatile liquid is spread on a solid surface; 
in technological applications, the liquid might be a paint or a lubricant. In this case, the liquid is neither in phase 
equilibrium with the solid substrate nor with the gaseous phase, and the total amount of liquid on the substrate 
is a conserved quantity 111. In the complete wetting situation, the liquid forms a continuous film on the substrate; 
in the partial or incomplete wetting case, the liquid forms droplets with the contact angle being determined by the 
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interfacial energies between the three phases meeting at the contact line |14|. Roughness of the substrate has been 
shown to cause contact angle hysteresis for advancing and receding contact lines |[L|Jlq| . 

A lot of work was specifically concerned with the interfacial and surface premelting properties of ice, due to 
its atmospheric and environmental consequences fl6| . Surface premelting of ice has been observed by a variety of 
experimental techniques ]T^-|l9|, and it is now believed that complete surface premelting (i.e., macroscopic growth 
of the surface liquid layer as the melting temperature is approached) only occurs for some orientations of the crystal 
surface and only if the vapor phase is diluted with air |2(J. Interfacial premelting of ice, giving rise to the low sliding 
friction of ice, has been deduced from wire regelation at low temperature [|2l|, and from viscosity measurements 
between surfaces of ice and quartz |2^] ; it also forms the basis for frost heave in frozen soils |23| . Complete interfacial 
premelting between ice crystals and a glass substrate was detected by ellipsometry The geometric structure 

of the glass surface was shown to play a vital role in this premelting phenomenon; in a series of experiments, the 
surface has been roughened by exposition to fluoric acid for different amounts of time, leading to surfaces with varying 
characteristic height-fluctuation amplitudes and wavelengths [^5| . For flat glass substrates the premelting was shown 
to be incomplete, while complete premelting was exhibited for glass substrates with a threshold amount of micro 
roughness |25| . 

The latter experimental observation motivated us to explore theoretically the possibility of a roughness-induced 
complete wetting or premelting transition. This describes the situation in which the flat substrate, for a given 
temperature, is not covered with a macroscopic liquid layer as coexistence is approached (corresponding to partial 
wetting), but, at the same temperature, is completely wet if the roughness of the substrate exceeds a certain threshold 
(in what follows we will use the wetting terminology both for the phenomena of premelting and wetting). In this 
paper we critically examine the conditions under which such a phenomenon can occur. As a result, roughness-induced 
wetting is indeed possible if the involved materials have the following properties: i) The tension of the substrate- vapor 
interface has to be larger than that of the substrate-liquid interface, ii) The surface area increase of the solid substrate 
due to its roughness has to exceed a certain threshold, which depends on the interfacial tensions of all three phases 
p5| . iii) The effective interaction between the two interfaces bounding the liquid layer for the flat case has to have 
a maximum for separations of about the amplitude of the substrate roughness. We also find that roughness-induced 
wetting is most likely to occur when the substrate roughness just exceeds a certain threshold value and will disappear 
for very large amplitudes of the roughness. 

The outline of this paper is as follows: In Sec. II we introduce the model and review some nomenclature for the 
case of a flat substrate. In Sec. Ill we extend the analysis to the case of a rough substrate. We first give necessary and 
rather general conditions under which roughness-induced wetting is possible. Using a formula which describes the van 
der Waals interaction between two sinusoidal surfaces, we then construct a sufficient condition for roughness-induced 
wetting in the form of a lower bound for the interaction between two flat interfaces. Finally, Sec. IV contains the 
discussion. 



II. FLAT SUBSTRATE 



Consider the situation as illustrated in Fig. 1, where a thin liquid film intrudes in between an inert solid surface and 
a top phase. The top phase can be either a vapor or a solid in thermodynamic equilibrium with the liquid film. The 
solid substrate, on the other hand, is completely frozen and far from its melting point. In this section, we will review 
results for an ideal solid substrate; namely, molecularly flat and homogeneous. Using the convention of labeling all 
physical quantities with a zero subscript for the flat case, the free energy per unit area can be written as 

F°(t)='YSL + 7 + P (£) + ^ (1) 

where 7s l denotes the solid-liquid interfacial tension, and 7 denotes the interfacial tension between the liquid and 
the top phase. The parameter /j, is the chemical potential difference between the liquid and top phases. Alternatively, 
it could also correspond to a Lagrange multiplier controlling the film thickness for non- volatile liquids with conserved 
total volume. The potential P°(£) represents the interaction per unit area between the two flat interfaces with a 
separation of £ and can be viewed as a thickness-dependent correction to the interfacial energies, depending both on 
the short and long-ranged parts of the molecular interactions [p6| . 

In the simplest approach, assuming pair-wise additive interactions between molecules and uniform densities in all 
coexisting phases, P°(£) can be expressed as 

/>oc p pO 

P°(£)= / dz d 2 p / dz' w(p,z-z') (2) 
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where w(p, z) corresponds to the local interaction energy difference per unit volume squared between the solid and 
the third phase. Dropping some constant terms, it can be written as 



w(r) = n 2 L U LL (y) - n L n s U LS (r) - n L n T U LT (r) + n T n s U T s(r) (3) 

Here, Uij(r) are the pair interactions between molecules and the n, are the particle number densities for each phase, 
where i and j are any of the relevant phases: solid (S), liquid (L), top phase (T) J5|. In a more realistic approach, 
one calculates (2) directly using the Lifshitz continuum theory of dispersion interactions |g7[ . 

In general, P° (£) is expected to vanish for £ — > oo as the two interfaces become decoupled and approaches a finite 
value for I — > . One therefore defines 

p°(i) = llt (4) 



v S for I -> 

where S is traditionally called the spreading coefficient and is given by 

S = 1ST - 1SL _ 7 (5) 

The above definition leads to .F (0) = 7st, just as one would expect: in the absence of any liquid, the total free 
energy is given by the interfacial tension between the solid and the top phase |30|. On the other hand, for an infinite 
layer of liquid coexisting with the top phase, one finds JF°(oo) = 75^ + 7 for /1 = 0, i.e., the two interfaces bounding 
the liquid do not interact and the total free energy is given by the sum of the two interfacial energies alone. One can 
notice that positive values of S = P°(0) — P°(oo) = JF°(0) — JF°(oo) correspond to a situation where an infinite liquid 
layer is energetically preferred over a vanishing liquid layer. Indeed, neglecting the possibility of additional minima 
of P (£) at intermediate values of £, positive and negative values of S correspond to wetting and non-wetting cases, 
respectively. On the other hand, in the case of a non-vanishing chemical potential [A, the minimum of the free energy 
will always be at finite film thickness, even for positive spreading coefficients' 

In the following, we will be exclusively concerned with the non-wetting case, i.e., S < 0. It will be convenient to 
modify the definition of the free energy slightly and to take the infinitely thick liquid layer as the reference state. The 
free energy difference, defined by A!F (£) = T°(£) — JF°(oo) and in the case of vanishing chemical potential, [i = 0, is 
given by 

Af°{£) = P°(£) (6) 

with the limiting values 

^ w = \lt a) 



v S for I -> 

Clearly, a wetting situation is realized if AjF°(£) > holds for all £ < 00 



III. ROUGH SUBSTRATE 

We introduce now the necessary framework to describe wetting on geometrically rough solids p^] . The free energy 
per unit projected area for a liquid film on a rough solid substrate (see Fig. 2) can be written as 

F{p, [Cl]) = v/l + [VC5(p)] 2 Isl + v'l+jvOW 7 + P(p, [a]) + MCl(p) - Cs(p)} (8) 

The solid and liquid surfaces are parameterized by Cs(p) an( l Cl(p), respectively, where p is a two-dimensional vector 
in a reference plane. The interaction term P(p, [Cl]) is a generalization of P°(£), Eq. (2), and is defined by 

poo r Ks(p + p') 

P(p,[Cl})= / dz / dV / dz' w(p',z-z') (9) 

JQl{P) J J-00 

Note that expression (9) is a local function of the liquid profile (l(p), but a non-local functional of the rough solid 
surface (s(p)- For the discussion of wetting behavior it is useful to average over the in-plane coordinate p, by which 
we obtain the effective wetting free energy 
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H£,[Cl]) = {Hp,Kl]))p (10) 

where we explicitly pulled out the dependence on the average film thickness 

i={C L (p)-c s {p)} p (ii) 

This parameter measures the average distance between the two interfaces. The effective free energy (10) can be 
written as 

T% [Cl]) = a slSL + a L (£, [Cl])7 + P(i, [Cl]) + ^ (12) 

In analogy to (10), the effective wetting potential P(£, [Cl]) is obtained from (9) by averaging over the in-plane 
coordinate 

p(e,[t L ]) = {p(p,[C L ]))p (i3) 

The parameters as and oli,{1, [Cl]) measure the ratios between the actual and projected areas of the substrate surface 
and liquid interface, respectively, and are defined by 



a s 



EE ( v/l + [VCs(p)] 2 )p (W) 



a l 



(£,[a])^(v/i + [va(p)p ) P (15) 



On the mean-field level considered in this paper, one can take the liquid interface to assume a fixed profile Cl(p) such 
as to minimize the free energy expression (12). By construction of the functional (12), this amounts to a constrained 
minimization of the free energy for a fixed average film thickness I. This yields the minimized free energy, denoted 
by T , as a function of £, 

T{t) = mm7{l^L])=7{l,[C L ]) (16) 
LCiJ 

The area ratio of the optimal liquid interface Cl has the limiting values 

al{l)^a L {t,[C L ]) = {T tl ( 17 ) 

since a very thin liquid layer follows the solid substrate roughness completely whereas a thick enough layer will be 
essentially flat (neglecting thermal capillary roughness). Just as for the flat case, P (£) = P(£, [££]) is expected 
to vanish for infinitely separated interfaces, i.e., P (oo) = 0. The interaction at contact (for £ — 0), is to a first 
approximation given by the interaction of the flat case times the surface area ratio of the rough solid surface, i.e., 
P*(0) re a s P°{0) = a s S ||. Defining the free energy difference by AT(£, [Cl]) = T(£, [Cl]) - 7*{oo), for which we 
set n = (the introduction of a non-zero chemical potential is straightforward and will be treated separately in Sec. 
III.C), one finds 

AT(£, [Cl]) = (a L (£, [Cl]) - lfr + Pft [Cl]) (18) 

The limiting values of the free energy (I) (obtained by minimizing with respect to the liquid interface profile Cl) 
are given by 

A?*m-/ Qs ^ ST_7sl )"^ for i ~ f0 (~\Q\ 
W - \ for £ -» oo [W) 

It is instructive to define the effective spreading coefficient, S e g = AJ 7 (0), which can be written as 

S c s = as(~fST - 7sl) - 7 (20) 
= (as - 1)7 + asS 

From the last equation it follows that the effective spreading coefficient is always larger than the bare spreading 
coefficient S, since 7 > and as > 1. The substrate area ratio can be expressed in terms of the spreading coefficients 
and the liquid interfacial tension as 

as = (21) 
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A. Definition of roughness-induced wetting 



With the definitions of the last sections we are now able to clearly define the subject and purpose of the present 
work. As already stated in the Introduction, we are concerned with the case where the flat substrate is not wet, i.e., 
the free energy difference AjF°(f) for the flat case, Eq. (6), is negative for some finite value of£. The central question 
is: Under which conditions will the rough substrate be wet, i.e., under which conditions does 

AT* (I) > (22) 

hold for all finite film thicknesses £ < oo? The answer of course imposes conditions both on the magnitude of substrate 
roughness (measured by as) and on the interactions between the coexisting phases, i.e., on the molecular interaction 
w(p, z), which enters the calculation of the wetting potential in (9). 

In Section III.B we give two rather general necessary conditions for roughness-induced wetting, which hold at very 
small wetting- layer thickness and rather large layer thickness (as compared to the substrate roughness amplitude), 
respectively. For the intermediate film thickness, we derive a sufficient condition in Section III.C. 



B. Necessary conditions for roughness-induced wetting 



1. Necessary condition for vanishing film thickness 



The necessary condition for a roughness-induced wetting transition which corresponds to (22) for vanishing film 
thickness (£ — > 0) follows from (19) and (20). It can be written as S e g > 0, which together with the non-wetting 
condition for the flat case (S < 0) leads to the inequalities j2j| 

1ST - 1SL < 1 < o>s(jst - Isl) (23) 

These inequalities can only be satisfied if -fsr > Jsl holds, since as > 1 by definition. From (5) one then also obtains 
that 7 > — S has to hold. Eq. (23) also imposes a lower bound on the substrate surface ratio as, 

as > TTsh (24) 

Clearly, the conditions (23) or (24), although necessary, are not sufficient for wetting, since the free energy AT {£) 
can develop a minimum at finite separation I. The conditions (23) and (24) simply correspond to a situation where 
the £ — ► oo thick liquid layer has a lower interfacial energy than a film of vanishing thickness I — ► 0. 



2. Harmonic approximation 

We now introduce the harmonic approximation, which we will adhere to in the remainder of this paper. Consider a 
corrugated solid surface, chosen to have sinusoidal undulations along one direction p\ of the two-dimensional reference 
plane (pi, p?) with amplitude hs and wave number q, 

C,s(p) = h s s\n(qpi) (25) 

The liquid profile is approximately (within linear response theory) characterized by the same q-mode undulation with 
a different amplitude hr,, vertically displaces by the film thickness £, 

( L (p) = h L sm(qp 1 )+£ (26) 

This geometry is depicted in Fig. 3. To linear order, there is also no phase shift between the two surfaces. In the 
following, the amplitude hs is assumed to be positive, with no loss of generality. The interfaces can not penetrate 
each other, constituting the non-crossing condition, which can be written as Cl(p) > (s{p), valid at any point p. 
This leads to the constraints 

h s ~£<h L <£ + h s (27) 

Expanding the expressions for the interfacial ratios as and a_L, (14) and (15) and keeping only terms up to quadratic 
order in the amplitudes and hs, leads to the following expressions for the area ratios 
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a s = l + \h 2 q 2 + 0((h s q) i ) (28) 



a L = 1 + -h\ q 2 + 0((h L qf) (29) 

which are expected to be good approximations for weakly corrugated interfaces (as long as h$q -C 1 and h^q -C 1). 
The necessary condition for roughness-induced wetting at vanishing film thickness (24) becomes 

as * 1 + \h 2 q 2 > Y^sff ( 3 °) 

We now identify two different physical regimes: for 7 ~ —S, defining the interaction dominated regime, the necessary 
condition for wetting leads to hsq 3> 1 and the expansion in terms of hq breaks down. Here, the solid roughness has 
to be quite pronounced and the behavior of the liquid interface turns out to be mostly dominated by the short-ranged 
part of the molecular interaction. For 7> -S, defining the tension dominated regime, condition (30) can be fulfilled 
even for small solid roughness (hsq <C 1); here, the liquid interface is dominated by its surface tension. It is the 
tension-dominated regime where the approximations leading to (28) and (29) and other simplifications made in the 
remainder this paper are valid; this is also the regime of most practical interest, since values for hsq characterizing 
rough surfaces in experiments are typically quite small |34| . 



3. Necessary condition for thick films 

An additional necessary condition for wetting is P (£) > valid for average film thicknesses approximately larger 

than the corrugation amplitude, £ ~ hs- This condition can be obtained in the following way: suppose we have a 
flat liquid interface, i.e., aj, = 1 • Then the first term in (18) vanishes, and in order for the minimized free energy 
difference AT (£) to be positive we have to require P (£) > 0. Clearly, a flat liquid interface is only possible for a 
liquid layer thickness which is larger than the amplitude of the solid roughness, otherwise the two interfaces sterically 
interact. Figure 4 schematically depicts the limiting case I « hs with the flat liquid interface just touching the solid 
substrate at the largest height fluctuation characterized by the amplitude /15. For sinusoidal interfaces described by 
(25) and (26) and using the non-crossing constraint (27), one obtains the inequality 

P*(£) > for h s < £ < 00 (31) 

For smaller distances, the interaction P (£) can actually become negative with AT (£) still being strictly positive, 
because then the free energy expression (18) always has a positive energy contribution from the interfacial tension of 
the liquid interface. 

This result has consequences for the important class of wetting potentials with a single minimum at finite but 
rather large wetting film thickness, which describe continuous wetting transitions as the minimum moves outwards 
to infinity: if the minimum occurs at distances larger than the roughness amplitude, it follows from (31) that the 
substrate roughness will not induce the wetting of the substrate. 

From the above considerations we see that the interaction P (£) has to have rather complex behavior; at zero £ or 
vanishing liquid film it is negative, since P (0) ~ asS and we start with the assumption of a non-wetting behavior 
(i.e., S < 0) for the flat solid surface. Only considering the £ = situation we see that there is a threshold value of 
the solid roughness in order to make the vanishing film limit energetically unfavorable compared to the infinite film 
limit, see (24). For a film thickness larger than the amplitude of the solid roughness, the interaction P (£) has to be 
repulsive in order to make roughness-induced wetting possible, see (31). The question that arises naturally is whether 
such a behavior is actually possible and what the conditions on the molecular interaction w(r) are. 



C. Sufficient condition for roughness-induced wetting 

In this section, we want to show for general film thicknesses under which conditions roughness-induced wetting, as 
defined by Eq. (22), occurs. In order to do so, we need to minimize the free energy with respect to the fluid interface 
profile Cl(p) for each average film thickness £ and for a given molecular interaction w(p,z) and interfacial tensions, 
according to (16) [with Cl(p) entering the expressions (8) and (9)]. We then have to check whether AT (£) > holds 
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for each £. If it turns out that AT* (£) < for a given £, we know that the system will prefer to have a stable film of 
this thickness, and wetting will not occur. 

For the present purpose, we turn this procedure somewhat around: we will try to determine the interaction for 
which roughness-induced wetting, for given values of the interfacial tensions and the substrate roughness, does occur. 
Furthermore, in many situations the molecular interaction w(p, z) is not easily available, and usually the wetting 
potential for two flat interfaces, P°(£), as defined by (2), is readily measured and calculated. We will therefore use 
a novel relation between the planar interaction P°(£) and the wetting potential between two sinusoidally modulated 
interfaces, P(p\,£, hi), which is derived in the appendix. Using this relation, we express the condition for roughness- 
induced wetting in terms of a lower bound on the planar interaction P {£). The analysis will be presented in the next 
two subsections. These parts are somewhat technical, and the unmotivated reader can easily skip these paragraphs 
and move on to the results in Section III.C.3. 



1. Construction of lower bound on P°(£) 

To proceed, consider first the range £ < hs- In the tension-dominated regime, defined by h$q <C 1, it follows that 
fg« 1 also holds. In this limit, the wetting potential, defined by (9), can for the special case of two sinusoidal 
interfaces (see Section III.B.2) be expressed in terms of the planar interaction P°(£). Neglecting curvature-like terms 
which turn out to scale like hsfiLq 4 , this relation is given by (see Appendix) 

p (w , W * (1 + w~W"/- < 32 » 

Clearly, the above form has the following desired property: for £ = one has hi = hs and one thus obtains for the 
spatially averaged potential P(£ = 0, hi) ~ asP°(£ = 0), as anticipated on intuitive grounds in the previous section. 
For either hs — or h^ = the formula (32) simplifies to P(p\,£, hi) = P°(£ + (hs — hi) sin[qpi]), which is exact. 

Still, the formula (32) is rather complicated and calculating the spatially averaged potential P(£,hr) = 
{P(pi, £, hi)) Pl is by no means trivial. The following observation is crucial: for a given value of t, the argument 
of P° in (32) is strictly smaller than 2£ (since 2hs£q 2 < 1). We can therefore choose a linear trial function of the form 

P$(t) = S + c(£)t (33) 

For a given value of £, we require the trial P^(t) to be a lower bound of the flat-interface potential P°(t), i.e., 
Pj-(t) < P a (t), in the finite range < t < 2£, which puts certain bounds on the slope c(£). It trivially follows that 
P(pi,£, hi) as given by (32) and evaluated with Pj. is always lower than P(p\,£, hi) evaluated with P° instead. 

The strategy will be as follows: If we can show that for a film with a given average thickness £ and using the 
linear trial function Pj, the resulting free energy difference AT (£) is positive, we know that the wetting film for this 
particular mean thickness is unstable. This follows since using P° instead of Pj. we will necessarily increase P (£) 
and thus also increase AT {£). If we succeed in showing the same for all values of £ (including £ = 0), it follows 
that P°(£) (which is by construction strictly larger than all the linear trial functions) is an interaction which shows 
roughness-induced wetting. In the subsequent calculations we actually revert this procedure and start with the linear 
trial functions, enforcing instability of the wetting film for each value of £ separately [leading to restrictions on the 
^-dependent slopes c(£)], and construct the function P°(£) from a superposition of all piece- wise linear functions. 
Choosing a linear trial function of the form (33) does not restrict the generality of our results, since one can express 
any function as the supremum of a set of suitable piecewise linear functions. 

Using the trial function (33) the averaging over the coordinate p\ can be easily done and,with (18) and (32), the 
resultant trial free energy is given by 

AT T {£, h L ) = -^hlq 2 + 5(1 + h s h L q 2 ) + c{£)£ + 0(q 4 ) (34) 
Minimizing this free energy expression with respect to hL, one obtains the following amplitude 

hi = -§^h s (35) 

which characterizes the minimizing liquid interface profile (~* L . The liquid interface corrugation is in phase with the 
substrate configuration, since S < 0. Since the two interfaces have to satisfy the non-crossing constraint (27), the 
minimizing amplitude h* L can be realized only for a thickness larger than some characteristic length £>£' . For £ = £', 
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the liquid interface has an amplitude given by (35) and touches the substrate surface; this situation is depicted in 
Fig. 5. This defines the length £', which is, combining (27) and (35), given by 

e' = h s (l + ^) (36) 

For smaller film thicknesses £, the amplitude Hl necessarily deviates from the value h* L and finally equals hs for I = 0. 
Note that the length £' is always bounded by hs/2 < £' < hs, as follows from the fact that 7 > — S, see (23). 



2. Calculation for S c s = 

In this paragraph we assume that S c s = 0, that means the substrate roughness just suffices to make the non- wetting 
state (I = 0) energetically unfavorable compared to the completely wet state (£ — 00) at coexistence (/j, = 0). This 
is obviously not necessary but renders the resultant expressions in a simpler form and is sufficient for fi = 0. The 
extension to general S e s will be done in Section III.C.4. From (30) the excess substrate surface area q 2 h 2 s /A is given 

by 



q 2 h 2 



s _ 



5 + 7 



For I < £' , the amplitude which minimizes the free energy and is in accord with the non-crossing constraint is given 
by hh = hs — £, that is, the interfaces just touch. Inserting this amplitude into the trial free energy expression (34), 
one obtains the minimal value (denoted by an asterisk) 



(2 



7 5 \ S£ + 



A ^ i) = c{i)e -^{s^) + Vs{-sTT) fOT e<e (37) 

Thus, the sufficient condition for the free energy difference to be positive is 

S fS + 2-f\ £ ( 

with the special values 



c(£) > -— — '- )+tt[ -r for £ < £' (38) 

h s \S + "fJ h 2 s\S + ^J 



m > -# (¥2) (») 



h s V S + 7 



c(0 > (40) 



and 

r(P'\ > - 

2h s V 5 + 7 

For £ > £', one inserts the expression (35) found for h* L into the free energy (34) and obtains 

AT* T {£) = c{£)£ + S- ^ 9 for £>£' (41) 

I67 



Using that 5 c ff = 0, it follows that 



Af* T (£) = c(£)£ + S+—— for £ > £' (42) 

4 7 (5 + 7) 

In this case, requiring the free energy to be positive leads to the condition 

^-35^ for l>e (43) 

with the limiting value (40) for £ = £' . 

Thus far, we have calculated a lower bound of the planar potential P°(t) for a given value of the average film 

thickness £, which consists of a linear function of finite extent (between and 2£) and ^-dependent slope. In order to 
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construct a lower bound on P°(t) which is valid for any film thickness, we have to calculate the upper envelope of 
these linear functions. This function we denote by P[° w (£). The sufficient condition for wetting, corresponding to the 
definition (22), can then be written as 

P°(£) > iCW (44) 

with Pi OVJ (£) being defined as 

P° w (£)=mzx{S + c(t)£} (45) 

The value of c(t) is given by (38) and (43) for t < £' and t > £', respectively. Since c{t\) > c(t 2 ) holds for any t\ < t 2l 
it is easy to see that the function Pj° w {£ ) can be expressed in closed form as 

P& w (£)=S + c(£/2)£ (46) 

Using the expressions (38) and (43), the function Pi OV! (£) is explicitly given by 

= {s(l- ggg) for 2£'<£< £ max (47) 



for e> 



max 



3. Results for S e g = 

Figure 6 shows the rescaled lower bound —P® ow /S as a function of the rescaled variable £/h$ for —S/j = 1/2. This 
function starts out at — f for £/hs — by definition of the spreading coefficient S. The slope at the origin is finite, and 
the function increases until it reaches its maximum at £ = 21' . We obtain a plateau in the interval 2£' < £ < £ max ; for 
£ > W the lower bound is given by — P® ow (£)/S = 0. In our estimate for the value of £ max we 

where £* is defined by P\ ow (£*) = (see Fig. 6). This is based on the observation that for £ > £ max the interaction 
as given by (32) is strictly positive as long as \Hl\ < h$, since in this case the arguments of P° in (32) are always 
larger than I*; consequently, the free energy difference (18) will be always positive. If, on the other hand, \hi\ > h$ 
holds, the free energy difference (18) will also be positive since then oil (£,h L ) > a L {0) and P {£) > P (0) holds for 
any £ > 0. 

The value for £* as obtained from (47) is given by 



£* S + 2 1 //5 + 2 7 \ 2 2(5 + 7) 



hs 7 V V 7 / 7 



(48) 



and reduces to the value I* /hs ~ 2 — \pl rj 0.59 in the tension-dominated regime, 7 3> — S, and approaches zero in 

the interaction dominated regime, 7 ~ —5. The largest possible value for l max /hs is thus l max /hs = 4 — \pl w 2.59. 
In the tension-dominated regime, 7 ^> — S, the assumption we started with, £q -C 1, thus holds for the whole range 
of thicknesses considered because the substrate roughness necessary to achieve roughness-induced wetting is small, 
hsq <C 1. In the same limit, the general expression (47) takes the values 



-S for 2<£/h s <4-\/2 (49) 

for £/h s > 4 - V2 



{ 

In the interaction dominated regime, for 7 ~ — S, the maximum of Pi ow (£) for 2£' < £ < £ max actually diverges as 
7 — ► —S, with 2£'/hs — > 1 and £ max /hs — > 2. Note that in this regime, also the substrate roughness necessary to 
induce wetting goes to infinity, sec (24). It follows that it is the tension-dominated regime where roughness- induced 
is most likely to be observed experimentally. 

The function P° ow {£) as given by (47) is plotted in Fig. 7 for -S/7 = 0, -S/j = 1/2 and -S/7 = 2/3, from bottom 
to top. As this ratio becomes larger, as one moves from the tension-dominated regime into the interaction-dominated 
regime, the maximum of —P® OVJ (i)/S increases, until it finally reaches infinity for —j/S= 1. 
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The interpretation of the results in Fig. 7 is the following: For a given ratio of the spreading coefficient S and the 
liquid interfacial tension 7, which are determined by the interfacial tensions of the problem alone and conspire to give 
S e e = 0, one rescales the amplitude of the wetting potential P (£) (which can be measured or calculated theoretically) 
by the spreading coefficient and the distance £ by the roughness amplitude hs- If the rescaled potential —P°(£/hs)/S 
is larger than —P® ow (£/hs)/S for all arguments, i.e. if the sufficient condition (44) is fulfilled, a complete wetting 
situation will result. 

An interesting feature of our results comes from the fact that we obtain universal lower bounds as a function of the 
film thickness I rescaled by the roughness amplitude hs- This suggests that there exists for a given substrate excess 
surface area ~ h s q 2 /A an optimal value of the roughness amplitude hs for roughness-induced wetting to occur. Since 
the maximum of the lower bound —P\ ow {£/hs)/ S is located at (./hs ~ 1, this optimal value of hs happens to coincide 
with the approximate location of the maximum of P°(£), the wetting potential. Assume that the sufficient condition 
(45) is in fact satisfied for this optimal value of hs- For roughness amplitudes much smaller than this optimal value 
(the function P° will be pushed to the right in the scaling plot Fig. 7) the sufficient condition cannot be satisfied for 
very small film thicknesses due to the finite slope of the wetting potential, suggesting an infinitesimally thin stable 
film (no wetting). For roughness amplitudes much larger than the optimal value (the function P° will be squeezed 
to the left in Fig. 7) the sufficient condition is not satisfied for average film thicknesses larger than some value at 
the order of the roughness amplitude, suggesting the formation of a very thin liquid film occupying preferentially the 
valleys of the substrate surface fluctuations (partial prewetting). 



4- Results for general S e s 

Here we treat the general case with S c s > 0, now including the cases where the substrate roughness as is larger 
than the necessary value obtained by setting S e s = in (21). In addition, we will require the free energy difference 
(18) to be larger than a constant denoted by Tq, as will turn out to be important if one looks at stable wetting layers 
in the presence of a chemical potential (see Section III.D). Using a trial function of the form (33) the following bounds 
for the slope c(£) follow 

effl > f {S + 2 Z } {SeS ~ S) -i " iS f ~ S) + for £<£' (50) 

h s S + f h% S + 7 £ 



Since c{£) for £ < £' now is non-monotonic and has a maximum at £ given by 



z =W (g 7f eff / o) ^ 

y 7(o e s - oj 

the global lower bound P® ow (£) can be calculated according to (46) only in the restricted range of 2£ < £ < 2£'. For 
£ < 2£, the function P® ow (£) as defined by (45) has a constant slope of c{£) given by (50) since £ < £' strictly holds. 
The global lower bound is thus given by 



S + h^§W) (( 5off - S K S + - 2 V7peff ~ FoKSeS - S)(S + 7 ) J for £<2£ 

I KTn Q a\ I 1 (■ s '+ 2 7)(Scff-5) I 2 7(Scff-S) f of < < Of 

,1 + Z{J-q ^cS)+ hs s+j fcf 2(S+ 7 ) ^S«S^« (53) 

S + 2^o + 5cffS 2 ^sS 27)2 for 2£'<£<£ max 

Fq for £ > £ max 



In Fig. 8 we plot P\ ow (£) for the fixed value —5/7 =1/2 and the values S c s = 0, S c «/S = — 1, and S e s/S = —2 (from 
bottom to top in the right portion in the graphs). The interesting feature is that the plateau of P® ov/ (£) for £ > 2£! 
increases with increasing S c s- This indicates that for a given wetting potential P (£) and a given roughness amplitude 
hs the roughness-induced wetting transition might disappear for wave numbers much larger than the threshold value 
determined by (30). 
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D. Stable films in the presence of a chemical potential 



So far, we have calculated a lower bound for the interaction, denoted by P^ w (£) , so that the free energy difference 
AlF*(£) is positive, Eq. (37), or larger than a positive constant T , Eq. (43), for all values of I. If P°(£) > P° ow (£) 
holds and the chemical potential vanishes, i.e., at coexistence, the substrate will be covered with an infinitely thick 
liquid layer. For non-zero chemical potential, however, the substrate will always be covered with a film of finite 
thickness. For clarity sake, let us assume the potential P°(£) to have a single maximum at finite separation, decay for 
large separations like P°(£) ~ a^ _CT (where a = m — 4 if the molecular interaction is given by (A7)) and be negative 
for zero separation, as corresponding to the non-wetting situation for a flat substrate. One notes that a = 2 for 
non-retarded van der Waals interactions. 

For finite chemical potential, the free energy then has a minimum at a finite separation £ m in- Assuming that 
£ m in ^S> q^ 1 , one can use Eq. (A21) and finds the leading terms of the free energy according to (18) 

AT(£, h L ) = IhW + P°(t) + \{h% + hl)P H {£) + (j1 (54) 

where the chemical potential term has been added. For a potential of the asymptotic form ~ a£~° this free energy is 
minimized by hi, — and the film thickness which is stable with respect to variations in £ is asymptotically given by 

£min ~ (aa/fi) 1 /^ (55) 

Choosing the constant Tq in (53) to be JF = AT (l m in) = AT{£ m i n ,hi = 0) and demanding P°(£) > P® ow (£) as 
given by (53), the minimum at £ m i n is indeed the global minimum and we have a stable wetting film of finite thickness. 
We will now do a local stability analysis for this free energy minimum. The linearized Euler equation, determining the 
stable liquid interface profile, can be obtained from the free energy expression (8) by requiring local force equilibrium, 
dF(p, [Cl]) I d(, L (p) = 0, and is given by 

7 ACl(p)+IIV)- Jd 2 p'{C L (p)~( S (p + p')-£}w(p',£) = Li (56) 

where U°(£) = — dP°(£)/d£ is the disjoining pressure for flat interfaces. Assuming the wetting layer to have a thickness 
corresponding to the global minimum as given by (55), which is equivalent to setting H°(£ m in) — A*, and approximating 
the interfacial profiles again by sinusoidal waves, Eqs. (25) and (26), the locally stable liquid interface amplitude turns 
out to be !| 

_ h s Jd 2 p cos(q Pl )w(p,£) 
L jq 2 +Jd 2 p w(p,£) 

Using formulas (A9) and (A25), this amplitude scale like h L ~ h s exp(-q£)£ 1/2 ^ m/2 q m/2 ^ 7 ^ 2 in the limit q£ > 1; 
for van der Waals forces (m=6) one obtains ~ hs exp(— q£)£^ 5 / 2 q^ x / 2 . In contrast to the global stability analysis 
leading to (55), where we averaged over the spatial coordinate p, the local equilibrium analysis now actually gives 
a non- vanishing amplitude h^. Inserting this amplitude back into the free energy expression (54), the stable film 
thickness is increased. However, the correction turns out to be less singular for the case of van der Waals forces and 
thus does not affect the asymptotic behavior of l m in as given by (55). 



IV. DISCUSSION 



Some necessary conditions for roughness-induced wetting have been obtained on very general grounds. First, the 
tension of the liquid interface has to be larger than the negative spreading coefficient, or, equivalently, the tension of 
the substrate- vapor interface has to be larger than the tension of the substrate-liquid interface. Second, the roughness 
of the substrate, measured by the excess area as compared to the flat substrate, has to exceed a certain threshold, see 
(24). Using an approach which is equivalent to a linear response analysis, we obtain a formula relating the molecular 
interaction between two rough interfaces to the interaction of two flat interfaces; the latter is the so-called wetting 
potential. Using this formula, a lower bound for wetting potentials in order to obtain this new type of wetting 
transition is derived. This lower bound constitutes a sufficient condition. More specifically, this lower bound has 
the following properties: i) for vanishing liquid film, or for zero separation between the two interfaces bounding the 
liquid layer, the interaction (which for this limit is the so-called spreading coefficient) is negative, corresponding to a 
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non-wetting situation in the case of flat interfaces; ii) the lower bound has a pronounced maximum with a height of 
at least the negative spreading coefficient at a separation of about the roughness amplitude. This height depends on 
the ratio between the liquid-interface tension 7 and the spreading coefficient S. 

These results concerning the possibility of a roughness-induced wetting transition agree with the very recent exper- 
imental observation of premelting induced by substrate roughness at the glass- ice interface |2q] , where glass surfaces 
showing a roughness with typical amplitudes of ~ Iran induced a complete interfacial premelting. Independently, 
the calculation of the dispersion interaction between half spaces of glass and ice separated by a liquid water layer 
show a pronounced maximum at about the same distance ~ lnm |p5| , consistent with our prediction for the lower 
bound of the wetting potential. We find the asymptotic behavior of the film thickness as a function of the chemical 
potential to be characterized by the standard van der Waals exponent, as shown in Sec. III.D. This is in disagreement 
with the experimentally measured exponent ~ 2 [p5[], which is to be contrasted with the exponent 1/3 as expected 
for non-retarded van der Waals interactions in the absence of any other interactions. Finally, we note that the in- 
terfacial energies for the case of interfacial premelting are influenced by grain-boundary energies in the disordered 
polycrystalline layers adjacent to the microscopically irregular substrate wall. In fact, it was argued that this 
effect could independently lead to a roughness-induced premelting transition. Such a mechanism could be described 
within our framework by assuming an interfacial energy 7 which depends explicitly on the roughness magnitude of 
the liquid interface. A further interesting effect might appear for the case of surface triple-point premelting: here a 
roughness-induced complete surface melting could be triggered by a roughening transition of the solid surface |?7| . 
Another interesting consequence of our results is that a discontinuous wetting transition might be converted into a 
continuous transition by a change of the effective wetting potential at small distances. 

Our calculations are most valid for small substrate roughness, as defined by hgq < 1; this also seems to be 
the experimentally most relevant limit, since rough surfaces produced by etching usually show modest corrugation 
amplitudes |34j . The closed- form expressions for the interaction of two corrugated interfaces obtained in the Appendix 
are to the best of our knowledge novel and applicable to a wide range of phenomena, including dewetting phenomena 
(where the whole argument has to be inverted, leading to roughness-induced dewetting). We also derive an expression 
for the curvature contribution to the free energy, which might play a role in the adsorption of membranes or vesicles 
on rough substrates. 
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The total interaction per unit (projected) area between the liquid interface at lateral position p located at Ck(p)> 
and a corrugated solid surface, parameterized by (s(p) (see Fig. 3), can be written as 



This expression is hard to deal with due to the non-local dependence on the shape of the solid surface, which enters 
in the integration boundary. Progress can be made by formally expanding (s(p + p') around (l(p) — The average 
separation between the two interfaces, I, is given by I = (Cl(p) — Cs(p)), where the brackets denote a spatial average 
over p. 

Keeping terms up to fourth order, one obtains 
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APPENDIX A: DERIVATION OF THE EFFECTIVE POTENTIAL 




(Al) 




dz / dV (Cs(p + p') - Cl(p) + i) w(p', z) + 
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I I A 2 p'^ s {p + p')-CL{p)+(.f§-w{p',z) 
l - J d 2 p>(( s (p + p')-( L (p)+t) 4 -^w(p',z) 



24 

where the expression for the interaction of two planar surfaces introduced in (2), 



(A2) 



P°(£) = f dzfd 2 pf dz' w(p,z-z') (A3) 

has been used. At this point it is useful to specify the interfacial profiles; we choose one-dimensional sinusoidal profiles 
for the liquid and the solid surfaces, as depicted in Fig. 3, as is sufficient and appropriate for a linear analysis, 

Cs(p) = Cs(Pi) = h s sm[q( Pl )] (A4) 

C L (p) = CL(pi) = h LS m[q Pl }+£ (A5) 

Now the expression (A2) can be averaged over the p-coordinates, thus yielding the mean interaction P(£, = 
{P(p, [Cl]))p, which for sinusoidal interfacial profiles reads 



P(£, h L ) =P\£) + -^J d 2 p (h 2 s +h\- 2h s h L cos[q Pl }) w(p, £) + 

— I d 2 p (h% + h\ + Ah 2 h 2 L - Ah s h L (h 2 + h\) cos[q Pl ] + 2h 2 h 2 L cos[2q Pl ]) g^w(p, z) 



(A6) 



and is an expansion up to fourth order in the interface modulation amplitudes hs and hL- 

To further proceed, it is appropriate to specify the molecular interaction w(p,z). In all what follows, an inverse 
power law defined by 

w(p,z)^A( P 2 + z 2 y m/2 (A7) 

will be used. Accordingly, non-retarded van der Waals interactions correspond to m = 6 with A being the Hamaker 
constant. The interaction between planar surfaces is given by 

P ° (e)= (m-2)(m-3)(m-4) ^ (A8) 
In addition, the following relations involving derivatives of P°(£) turn out to be useful: 

p(U)(l) = ^_ £ 2-ra = A fd 2 p (p 2 +£ 2 )- m / 2 (A9) 

m - 2 J 

P {IV \£) = 2nA(m - 1) £~ m (A10) 

y > (m-2)(m-3)(m-4)(m-5)(m-6) 1 ' 

where d 2 P { - H \£)/d£ 2 = P°(£). Note that P^ n \£) is not defined for van der Waals interactions with m = 6; this 
important case will be considered separately. 

The terms in (A6) involving a cosine can now be evaluated analytically; for the first term, and using the interaction 
defined in (A7), we obtain 

[°° [°° cos[g Pl ] AV^T(^) f°° cos[q Pl ] 

A / d ?2— — = 7T7^ / dPl 



\pl+pl+( 2 ) m ' 2 ~ ^f) J-oo ' (pl+pf^' 2 

= T^^-iMX^r 72 ^ 2 ^ 1 -" 1 / 2 (A12) 
In the last equation, K v (z) denotes the Modified Bessel Function in standard notation. 
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1. Small separations 



Let us first concentrate on separations i much smaller than the typical wavelength q 1 of the surface corrugation, 
i.e., I -C q ■ For small values of the argument z, the product K v [z)z v can be expanded as a power series 



K v {z)z v — ao + a,2Z 2 + a4Z i + ■ 
with v = m/2 — 1 and the coefficients given by 



a = 2^ 1 I» 



for v > 



a 2 



of— 3 



2"- J r(^ - 1) for 1/ > 1 



a 4 = TST(v - 2) for i/ > 2 
Note that for van der Waals forces, characterized by v = 2, the coefficient 04 is given by 

a 4 = (3/2-2 7 + 21og2-21ogz)/16 



(A13) 



(A14) 
(A15) 
(A16) 



(A17) 



with 7 being the Euler constant defined by 7 = 0.57721. Using this expansion and the definitions (A8-A11), the 
cosine term in (A6) can be written as 



.4 



dpi 



dp2 



cos[qpi 



P (")(|) 



m "3j D o^ , ( m ~ 3)(m- 5) 4 ( _ Jf) 



For the case of van der Waals interactions, m = 6, one analogously obtains 



g*pl-"J(*) + 0(g"') (A18) 



dpi 



cos[gpi] 



2tt ( 1 



£ 4 4£ 2 



+ 



3_4 7 + 41n2-41n(^) 
64 ' 



0(<Z 6 ) 



(A19) 



which includes a logarithmic singularity for small separations. 

For the other term in (A6) involving a cosine one can interchange the differentiation and integration (for m > 6); 
the additional integral needed involves cos[2qpi] and is given by 



A 



dpi 



cos[2gpi] 
dP2 7^ 



P {II \l) - 2(m - 3)q 2 P°(£) + 2(m - 3)(m - 5)g 4 pf-^(^) 



(A20) 



Using the formulas (A18), (A20), and the definitions (A8)-(A11), one obtains the following expansion for (A6) 

(m — 3) (to — 5) 



P(£,h L ) = P°(£) 

+P (II) {t) 



m — 3 , , 9 
1 + ——h s h L q 2 



128 



-^(/,5-^)V + ^ m " 3)(m " 5) ^^" 4 



G4 



-h%hitt 



+ P {IV) {^{hs - h L f - 5P(t, h L ) + 0(h 6 , q 6 ) (A21) 
64 



which is valid for q£ <^ 1 and to > 4 (for to = 4 additional logarithmic singularities appear in terms proportional to 
q ) , The terms up to 0(h 2 , q 2 ) have been obtained previously in the context of the dynamics of thin wetting layers 
p8[. The correction term SP corresponds to a curvature contribution and is given up to 0(h 4 , q A ) by 



sp(i,h L )= {m - 3 ] ( r- 5) p^ i Ht)h S h^ 

16 



(A22) 



This curvature contribution has additional non-analytic terms for the case of non-retarded van der Waals interactions, 
to = 6; for this case, one obtains 



KT>(P h \ n 3-4 7 + 41n2-41n(gl) L 4 
dP(£,h L ) = — h s h L q 



256 



(>4 



(A23) 



The same singularity ~ q 4 In q has been found for a free interface in the presence of van der Waals forces. E6 
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2. Large separations 



In the other limit, for q£ » 1, the product K 1/ {z)z v is given by 

K v {z)z v = r^e-^l + (4^ 2 - l)/8z + • • •) (A24) 

In this case, the expression (A12) is asymptotically given by 

/°° r°° rr.i\nn.] o3/2-to/2_3/2 a 

dpi / d P2 gggbg = 2 r * ^ ^-./^i/^/V (A25) 

From equation (A6) one immediately obtains 

P(4 fc L ) = PV) + Jp (/7) W(4 + *£) + ^P (IV) W(hs + hi + 4h%h 2 L ) + 0(e-^) (A26) 
which is valid for » 1. 

3. Closed-form expressions for the effective potential 

In the following, formulas are presented, which express the series for the effective interaction between the two rough 
interfaces, (A21) and (A26), in terms of the interaction P°(£) between two flat interfaces. 
For the case ?f<l, this expression is given by 

P(pi,l,h L ) = (l + ^WcosVi]) 1/2 ^° ( f^rf ~ 2 kL) P^rk ) SPfal, h L ) (A27) 

For the case ql 3> 1, the corresponding expression is given by 

P{p u £, h L ) = (P° {£ + h L sinfepi] + h s sin[gri])) ri + 0{e~ qe ) (A28) 

where T\ is the local lateral coordinate on the liquid interface and is averaged over, leaving only the dependence on 
the coordinate p\ in the substrate interface. Expressions (A27) and (A28) depend explicitly on the spatial coordinate 
Pi, that they indeed reproduce term by term the series (A21) and (A26) can be checked by expansion and averaging 
over pi . The validity of the closed form expression is thus proven for power laws with arbitrary m; we were not able to 
extend this proof to interactions which include a cut-off at small separations. However, it is likely that the formulae 
(A27) and (A28) are also accurate for potentials P°(£) which do not diverge as £ — > 0. This is supported by the fact 
that for £ = the formula (A27) exactly describes the surface-like energy contributions. 
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FIG. 1. Schematic view of a flat substrate: A liquid layer of thickness £ intrudes between the inert substrate and the top 
phase, which can be either the vapor or the solid in chemical equilibrium with the liquid, corresponding to wetting or interfacial 
premelting, respectively. The liquid-substrate and liquid-top phase interfacial energies are denoted by jsl and 7, respectively. 

FIG. 2. Liquid film on a rough substrate: For thin mean film thickness £, defined by the averaged local separation between 
the two interfaces, the liquid interface follows the substrate corrugations. 

FIG. 3. Simplified geometry in the single q-mode approximation, with the substrate surface parameterized by Cs(pi), an d 
the liquid interface parameterized by Cl(pi), shown along the direction parallel to the wavevector of the sinusoidal profile. The 
mean separation I corresponds to the distance between the mean positions of the interfaces, denoted by broken lines. 
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FIG. 4. Minimal mean film thickness I for which a flat liquid interface is still possible. The liquid interface touches the 
substrate surface at isolated points, and the film thickness i corresponds to the characteristic corrugation amplitude hs- 

FIG. 5. Liquid film with a thickness £' , for which the liquid interface with a corrugation amplitude h* L as given by (35) just 
touches the substrate. Since the two interfaces cannot cross, in the single q-mode approximation the amplitude hi, increases for 
smaller values of I until one finally obtains hi — hs in the limit I — > 0. In a more realistic model, one might obtain a ruptured 
interface for small layer thicknesses, corresponding to separate droplets. 

FIG. 6. Plot of the lower bound for wetting potentials, —P^/S, as a function of the rescaled film thickness £/hs, for 
— S/7 = 1/2 and S e g = 0. In the interval < t < 21', P\ OVJ {i) is a monotonically increasing function. In the interval 
21' < I < < m », -P/owW is constant. For I > ^ max one finds P° ow (£) = 0. The value t is defined by P° ovr (t) = 0. 



FIG. 7. Plot of the lower bound for wetting potentials, —Pum/S, as a function of the rescaled film thickness £/hs, in the 
limit 7 3> — S, for —S/j = 1/2, and for —S/j = 2/3 (from bottom to top). The curves shown correspond to the special case 
S e s = 0, i.e., the relative area of the substrate is just sufficient to make the free energy of the dry state (£ = 0) higher than the 
completely wet state [I = oo). The area ratios are given by as = 1, 2, and 3, from bottom to top. 

FIG. 8. Plot of the lower bound for wetting potentials, —P^/S, as a function of the rescaled film thickness l/hs, for the 
fixed value —S/j =1/2 and for S e g = 0, S c e = — S, and for S e s = —2S (from bottom to top for the right portion of the plots). 
The relative area increase of the rough substrate is given by as = 2, 3, and 4, from bottom to top. 
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Fig.4, Netz and Andelman 
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Fig. 5 , Netz and Andelman 



Fig. 6 , Netz and Andelman 
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Fig. 7 , Netz and Andelman 
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Fig. 8 , Netz and Andelman 
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